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CHAPTER I 


INTRODUCTION 

The studies reported here were part of a long-range program conducted 
in the UCLA Molecular-Beam Laboratory with the objective of adding to the 
available knowledge concerning energy and momentum transfers between 
satellites and upper-atmosphere gases. Important applications of such 
knowledge include (a) estimating the lifetimes of satellites for given 
atmospheric models and (b) extracting upper-atmosphere densities from 
satellite-drag data. 

The aerodynamic drag of a satellite is predicted sometimes using a 
model in which the energy accommodation between the incident molecules and 
the satellite surface is complete and the reflected molecules are scattered 
diffusely (with a cosine scattering distribution). However, deviations 
from this idealized model could yield significantly different satellite 
drags. Hence, specific objectives of the present studies included (a) 
measurements of energy accommodation as a function of incidence angle and 
scattering direction, (b) measurements of scattering distribution as a 
function of incidence angle, and (c) calculation of the drag coefficient 
from the forementioned measurements. 

The studies of the interactions of satellite-speed helium atoms with 
satellite-type surfaces were initiated by Dr. S. M. Liu. Spatial distri- 
butions of helium atoms reflected from several satellite-type surfaces 
are summarized in Ref. 1; energy distributions of helium atoms reflected 
from a 6061 T-6 aluminum surface are summarized in Ref. 2. In the present 
report, spatial and energy distributions of helium atoms scattered fx ,m an 
anodized 1235-0 aluminum surface (also tangential and normal momentum 
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accommodation coefficients calculated from these distributions) are reported 
Included are (a) a procedure for calculating drag coefficients from measured 
values of spatial and energy distributions and (b) the drag coefficient 
calculated for a 6061 T-6 aluminum sphere from the data reported in Refs. 1 
and 2. The drag coefficient for an anodized aluminum sphere will be 
reported in a later publication. 
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CHAPTER II 


EXPERIMENTAL METHODS 

The present study was carried out in the UCLA Molecular-Beam Labora- 
tory using the molecular-beam system shown schematically in Fig. 1. The 
satellite-speed helium beam was generated using the arc-heated supersonic 
beam source developed by Young ^ and passing the central portion of the 
arc-heated jet through a skimmer and two collimating orifices. The re- 
sulting beam was characterized by a multi-disk velocity selector whose 
orientation and angular speed were adjusted to yield a 7000 m/sec beam. 

The detection system (Fig- 1) facilitated measurements of the three- 
dimensional scattering distributions and the energy distributions of the 
scattered atoms. The system included a target positioning mechanism, a 
detector rotating mechanism and a mass spectrometer and/or a retarding- 
field energy analyzer. These components are described in Ref. 2. 

Spatial -Distribution Measurements 

Spatial distributions of a helium beam reflected from an anodized 
aluminum surface were measured essentially as described^in Ref. 2. The 
main difference in the present measurements was that the distance between 
the surface and the entrance orifice of the ioniser was reduced to one 
inch (from the two inches used previously) in order to increase the signal 
intensity. A block diagram of the electronic system for the mass spectro- 
meter and the multi-disk velocity selector is shown in Fig. 2. 
Energy-Distribution Measurements 

The energy distributions were measured using the retarding-field 
energy analyzer described in Ref. 2 and the detection system (a slight 
modification of the one described in Ref. 2) shown schematically in 
Fig. 3. The energy spectrum of the reflected atoms at a given scattering 




Figure 1 , Schematic Diagram of the Molecular Beam System. 














quadrupole surface 
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Figure 2. Block Diagram of the Electronic System of the Mass Spectrometer 
and the Multi-Disk Velocity Selector. 




































angle was obtained by measuring the reflected-beam density as a function 
of the retarding potential. In the arrangement described in Ref. 2, the 
incident beam was not modulated so that the measurement included contri- 
butions from both the reflected beam and the background; the background 
contribution was determined in a separate measurement by intercepting the 
incident beam by a flag located inside the detection chamber. In the 
arrangement used here, (a) the beam entering the detection chamber is 
modulated by a chopper rotating at 10 cps before it impinges on the target 
surface, (b) the signal from the SSR counter is fed to an Ithaco Dynatrac 
391A lock-in amplifier and (c) the output of the lock-in amplifier 
represents the contribution from the beam alone since the background contri-^ 
bution is subtracted automatically. 

Surface Preparation 

The target used in the present studies was cut from a 0.005-in. -thick 
anodized 1235-0 aluminum foil provided by personnel at the NASA Langley 
Research Center. The surface was cleaned with acetone prior to installa- 
tion on the target holder. 
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CHAPTER III 


EXPERIMENTAL RESULTS 

Spatial Distributions 

Spatial distributions measured for a satellite-speed (7000 m/sec) 
helium beam scattered from an anodized 1235-0 aluminum surface for six 
different incidence angles (0°, 15°, 30°, 45°, 60 °, and 75° from the sur- 
face normal) are shown in Figures 4-9. In these figures, the center of 
the polar diagram corresponds to the point of impingement. The incident 
beam impinges on the test surface (which coincides with the surface of 
the page) from the bottom of the diagram with the given incidence angle 
measured from the surface normal. The upper (0^ > 0) and lower (0 < 0) 

halves of the diagram represent the forward-scattering and backward- 
scattering regions respectively. The dashed lines at constant value of 
0 r indicate detector paths from <j> = 0° to <j> = 90°, where <j> denotes the 
out-of-plane scattering angle. 

These plots exhibit trends similar to those found for a 6061 T-6 
aluminum surface [2], Of special significance (for both the present 
surface and the 6061 T-6 surface) is the backscattering, which becomes 
more prominent as the incidence angle increases toward the surface 
tangent (i.e., for large values of 0^). Such backscattering yields relati- 
vely high drag coefficients. 

E nergy-Accommodation Coefficients \ 

Measurements of energy distributions of s&tellite-speed helium atoms 
scattered from an anodized aluminum surface were made for six different 
incidence angles (0^ = 0°, 15°, 30°, 45°, 60° and 75° from the surface 
normal) . For each incidence angle, distributions were measured at approxi- 
mately forty scattering positions. These scattering positions were chosen 
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from eleven in-plane scattering angles (0 = ± 75 e , ± 60°, ± 45°, ± 30°, 

± 15 c and 0°) and six out-of-plane scattering angles ( 4 . = 0°, 15°, 30°, 
45°, 60° and 75°). Measurements were not possible within a solid angle 
around the incident beam (due to interference between the detector and 
the incident beam at these scattering positions) and for some glancing 
scattering angles (due to weak signal-to-noise ratios). 

The data obtained for the net helium beam as well as for the back- 
ground helium was recorded both on chart-recorder paper as well as on 
computer cards. It was analyzed using the procedure described in Ref. 2. 
The resulting differential energy accommodation coefficients (AC) E ( 6 ^, 

6 , <£0 obtained at all possible scattering angles are given in Tables 
I -VI. 

The differential energy accommodation coefficients obtained show some 
fluctuations, due apparently to the weak signal-to-noise ratio which 
results from the relatively diffuse scattering. 
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TABLE -1 . THE p»f FE.E6NTIAL ENERGY ACCOMMODATION 

COEFFICIENTS FOR 1 OOO m /sec HELIUM 
BEAM SCATTERED FROM A^ODUED ALUMINUM 
SURFACE AT 0° INCIDENCE ANGLE. 
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TABLE -UL . THE DIFFERENTIAL ENERGY ACCOMMODATION 
.COEFF/CIEWTS’ FOR 7 OOO m jsec HELIUM 
BEAM SCATTERED FROM AUODHED ALUMINUM 
SURFACE AT 30° INCIDENCE ANGLE. 
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TABLE -3T . THE DIFFERENTIAL ENERGY ACCOMMODATION 
COEFFICIENTS FOR 7000 rn /see HELIUM 
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TABLE -"2C. THE DIFFERENTIAL ENERGY ACCOMMODATION 
.COEFFICIENTS' FOR 700O rn /sec HELIUM 
beam scattered from anodized aluminum 
SURFACE AT 15 * INCIDENCE ANGLE. 





























CHAPTER IV 


EVALUATION OF DRAG COEFFICIENT FOR THE 6061-T6 ALUMINUM SURFACE 

This chapter summarizes the procedures used to deduce a drag co- 
efficient for spheres from measured differential energy^accommodation 
coefficients, and gives the value deduced for a 6061 -T6 aluminum sphere 

( 2 ") 

using the measurements made by Liu 

Extrapolation of Scattering-Distribution and Energy-Accommodation 
Measurements 

Recall that scattering and energy distributions could not be measured 
in a significant portion of the back-scattering region due to interference 
of the detector with the beam. In addition, energy distributions could 
not be measured for large values of the scattering angle, 6^, and the 
azimuthal angle, 4>> due to the weakness of the signal at these angles. 
Hence the following method was devised to extrapolate the measurements 
to regions where measured values are not available. 

If, for a given value of 0. , one plots the scattering distributions 
measured in two orthogonal planes (one plane defined by a constant value 

i : 

of <p, the other by a constant value of 0 ) , then the scattered-beam 
density at the intersection of these two planes must be common to both 
planes. Thus, the density at a particular point could be determined, if 
the values in its vicinity in two orthogonal planes passing through this 
point were known, by extrapolating in both planes to a common value at 
this point. The estimated value at this point could be used now in the 
estimation of the density at a new point in its vicinity. 

This procedure was used to fill in the gaps in Tables III-l - II I -6 
of Ref. 2, where measured values of density and energy accommodation 
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coefficients are reported. Both the measured and extrapolated res;. ,s are 
given in Tables VII - XII of the present report. 

Differential Normal and Tangential Momentum Accommodation Coefficients 
Differential momentum accommodation coefficients may be defined as 


(Pi * PJ 


r'N 
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(Pi - P r ) T 


tv l - Vt 
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where a NM and a^ M are respectively normal and tangential differential 
momentum accommodation coefficients, subscript N refers to the normal 
component and subscript T to the tangential component. The velocity 
components are given (cf. Fig. 10) by 
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The fraction of particles leaving in solid angle do (= u sin0d<j>du) is 
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TABLE - JX . 


THE DIFFERENTIAL ENERGY ACCOMMODATION COEFFICIENTS AWL 
THE SPATIAL DENS IT y DISTRIBUTION* FOR 7000 rnjsec HELIUM 
BEAM SCATTERED FROM CLEANED 6061 - T6 ALUMINUM 
PLATE AT 30 ° INCIDENCE ANGLE . 



(a) DIFFERENTIAL ACCOMMODATION COEFFICIENT (Y) 

(b; SPATIAL DENSITY (ARBITRARY UNITS) 

^MEASURED VALUES ARE UNDERLINED. 
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(a) DIFFERENTIAL ACCOMMODATION COEFFICIENT (X) 
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The direction of v also is opposite to that of v. . 

W N X N 

The computed values of the differential normal and tangential 
momentum accommodation coefficients are given in Tables XIII - XVIII. 
Overall Normal and Tangential Momentum Accommodation Coefficients 

Overall normal and tangential momentum-accommodation coefficients 
a were calculated for each incidence angle by averaging over the differen- 
tial coefficients cu taking into account the flux distribution and the 
solid angle represented by each measurement. Thus, the overall accommo- 
dation coefficient was deduced from the differential coefficients using 


a = 


■?j . a. (Aw) . 
O . (Aw) . 

i J i 


where is the flux and (Aw)^ is the solid angle corresponding to the 


measurement. Since is proportional to n^ v^l-a^ and (Aw)^ to cos^. 


one may write 

£ n.Jl-ac a A cosiji. 

i c . 1 1 

— l l 

a = 


. n. Jl-a« cosd>. 
l E. i 

l 


where n^ is the number density, is the energy-accommodation coefficient, 
and is the out-of-plane scattering angle measured from the plane of 
incidence. 

The deduced values are plotted as a function of incidence angle 9^ 
in Fig, 11. 
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TABLES XIII 


XVIII 


DIFFERENTIAL NORMAL AND TANGENTIAL MOMENTUM 
ACCOMMODATION COEFFICIENTS 


In the following tables: 

otg = differential energy accommodation 
coefficient 

n = spatial density 

= differential normal momentum accommoda 
tion coefficient 

= differential tangential momentum accom 
modation coefficient 
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Evalu ation of Drag Coefficient fo r Spheres 

An expression for the drag coefficient for spheres moving in a rare- 
fied atmosphere was derived in terms of normal tangential momentum- 
accommodation coefficients. (See Appendix A.) For the case in which the 
molecular random thermal speeds in the atmosphere are negligible in com- 
parison with the satellite speed, one obtains 


4 J / 2 r 

'D = CV V w > N (a NM c0s9 i ) + (V c O -V T ra TM Sine i ) 

i 


x sin6. cos 9 . A8. 
1 li 


where and are, respectively, overall normal and tangential momentum- 
accommodation coefficients (a function of 6.), V is the incident speed of 

2. oo * 

the beam for which the experimental measurements are made, and is the 

average reflected speed for complete accommodation with tne surface. 

The drag coefficient for 6061-T6 aluminum spheres was computed from 

the above expression using the results in Fig. 11 by numerical integration. 

(A computer program for this calculation is given in Appendix B.) For 

the given conditions (He impinging at 7000 m/sec), = 2.b4 was obtained. 

The derivation of an expression that takes into account the mclecular 

random thermal speeds in the atmosphere is more difficult. However, for 

the particular case when the normal and tangential momentum-accommodat ion 

coefficients are constant, a relatively simple expression for the drag 

coefficient may be derived in terms of the speed ratio, S. See, e.g., 

(4) 

p. 167, Equation (16) of Patterson . This predicted dependence on S, 

and some experimental points obtained for He, are shown in Fig. 12 (from 
Ref. 4). The good agreement between the predicted and the measured values 
establishes the trend of the variation of the drag coefficient with the 
speed ratio. Note also that in Fig. 12, for S greater than about 2.4, 


40 



SPEED RATIO, S 


FIG. 12. Drag Coefficient as Function of Speed 
Ratio (from Fig. 38 of Ref. 4). 
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the drag coefficient approaches a terminal value and becomes relatively 
independent of the speed ratio. 

For an atmospheric temperature of 1000°K, the c ^eed ratio for a 7000 
m/sec incident helium beam is 3.44. Thus, the expression for the drag 
coefficient given here is sufficiently accurate for many satellite 
conditions . 
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APPENDIX A 


Deriva tio n of Drag Coefficient f or Spheres for the Case of 
Inf in j.teSp ee d Ratio 

Consider a sphere of radius r which intercepts a beam of molecules 
with uniform speed V as shown in Fig. 15. Expressing the force exerted 
on a toroidal differential area dA as a function of the average normal 
and tangential momentum accommodation coefficients (functions of 9.) and 
summing over the upper half area of the sphere, one obtains 


V no • 0 f c olli s i ons 
time x area 


dA 


( V I ’wV c ‘NM c 0 se i ) * <P i -Pw ) T (a TM sin V 


1 2 
= t C. nm V A 
2D ® 


The number of collisions per unit area per unit time is nV^; the differential 
area is 

2 

dA = 2nrsin9 . (rd6 . ) z os6. = 2irr~ sine. cos9.d9. 

ill ill 

Then 


X nV 2zt~ sin9.de . 

^ ® li 


'Pi-PwVNM^V^Pi-PwVW^V 


cos9 . 

l 


= -r C. nm V “A 
2D °° 


which may be simplified to 


!tt V r~y 


(P 


• -P ) » (*vm cos9 ; ) + (P- -p )_(ot_.,sin9. ') sin9 . cos© . d9 . 
l *VN NM i *i r w'T TM i J i li 


“ 2 C D m V A 

2 

Since A = rr~ and p = mV, the drag coefficient C^ may be written 


r = 
"D 


4 */ 2 r 1 

= — y (V -V ) „a.,,.cose . + (V -V 1 a sin9. x sin© . cos© . d9 . 

V 1 » w'N NM i “ w T TM l l l l 

00 9 . =0 L -I 
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Fl£. 13 REPRESENTATION OF A SPHERE MOVING IN A 
RAREFIED ATMOSPHERE. 


45 


For the co-ordinate system shown in Fig. 10, 
V = V cose. 

oo ao 1 

vi 


V = V sine. 

ao oo 



V = o 

W,. 


The latter five equations were used in the evaluation of the drag co- 
efficient given in Chapter IV of this report. 
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